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The first observation of the production of a W boson with a single charm quark (c) jet in pp 
collisions at y/s = 1.96 TeV is reported. The analysis uses data corresponding to 4.3 fb _1 , recorded 
with the CDF II detector at the Fermilab Tevatron. Charm quark candidates are selected through 
the identification of an electron or muon from charm-hadron semileptonic decay within a hadronic 
jet, and a Wc signal is observed with a significance of 5.7 standard deviations. The production 
cross section awe (pre > 20 GeV/c, |7j c | < 1.5) x B(W — > iv) is measured to be I3.6I31 an d 
is in agreement with theoretical expectations. From this result the magnitude of the quark-mixing 
matrix element V cs is derived, \V CS \ = 1.08 ±0.16 along with a lower limit of |V CS | > 0.71 at the 95% 
confidence level, assuming that the Wc production through c to s quark coupling is dominant. 
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The associated production of the W boson with a sin- 
gle charm quark in proton-antiproton collisions is de- 
scribed at lowest order in the standard model (SM) by 
quark-gluon fusion (gq — > Wc), where q denotes ad, s, 
or b quark. At the Tevatron proton-antiproton collider, 
the larger d quark parton distribution function (PDF) 
in the proton is compensated by the small quark-mixing 
(Cabibbo-Kobayashi-Maskawa or CKM) matrix element 
\V c d\, so that only about 20% of the total Wc production 
rate is due to gd — ¥ Wc, with the majority due to strange 
quark-gluon fusion. The contribution from gb —> Wc is 
also heavily suppressed by \V c b\ and the b quark PDF. 
The Wc production cross section is therefore particularly 
sensitive to the gluon and s quark PDFs P, Q , at a mo- 
mentum transfer Q 2 of the order of the W boson mass 
(Mw), and to the magnitude of the CKM matrix element 
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V cs . Measurements of Wc production in high energy pp 
collisions are of interest because they constrain the pro- 
ton's s quark PDF at momentum transfers about three 
orders of magnitude higher than in neutrino-nucleon scat- 
tering Q. Finally, the Wc final state is similar to final 
state of other processes, such as single top-quark produc- 
tion, neutral and charged Higgs boson production, and 
supersymmetric top-quark production. The techniques 
developed here could lead to a better understanding of 
those samples and their searches. Calculations of W + 
heavy quark production are available at leading order 
(LO) and next-to-leading order (NLO) in quantum chro- 
modynamics (QCD) [4], with the NLO cross section pre- 
diction about 50% larger than the LO calculation. Over- 
all, the uncertainty on the NLO theoretical expectation 
for the Wc production cross section at the Tevatron is 
10-20%, depending on the charm phase space considered, 
dominated by uncertainties from the choice of factoriza- 
tion and renormalization scales and the shape of the s 
quark PDF. 

We present the first observation of pp — > Wc produc- 
tion. The charm quark is identified through the semilep- 
tonic decay of the charm hadron into an electron or muon 
(referred to in this Letter as "soft leptons"). This mea- 
surement supersedes our previous result [5[, where the 
cross section for pp — > Wc was determined with a preci- 
sion of approvimately 30% and a statistical significance 
of about 3 standard deviations. The present analysis 
is performed using a data set more than twice as large 
and signal events with soft electrons are included to in- 
crease the acceptance. The analysis exploits the correla- 
tion between the charge of the W boson and the charge of 
the soft lepton from the semileptonic decay of the charm 
hadron. Charge conservation in the process gq — ¥ Wc 
(q = d, s) allows only W + c and W~c final states; as a re- 
sult the charge of the lepton from the semileptonic decay 
of the c quark and the charge of the W boson are always 
of opposite sign, neglecting any effects due to slow-rate 
charm quark oscillations. 

The W boson is identified through its leptonic decay 
by looking for an isolated electron (muon) carrying large 
transverse energy Et (momentum pr), with respect to 
the beam line. The neutrino escapes the detector, caus- 
ing an imbalance of total transverse energy, referred to 
as "missing Et' ' ($t) dl- Quarks hadronize and are ob- 
served as jets of charged and neutral particles. Charm 
jets are identified by requiring an electron or muon can- 
didate within the jet ("soft lepton tagging" or "SLTf"). 
Events are classified based on whether the charge of the 
lepton from the W boson and the charge of the soft lep- 
ton are of opposite sign (OS) or same sign (SS). The 
Wc production cross section is then calculated using the 
formula 



where iV^f ~ ss {N°^~ ss ) is the difference in the num- 
ber of OS and SS events in data (background), A is the 
product of the efficiency, for identifying Wc events, with 
the kinematical and geometrical acceptance, and J L dt 
is the integrated luminosity of the data sample. The 
quantity S = (N° s c - N^ S C )/(N^ S C + N SSj accounts for 
the charge asymmetry of the sample of real reconstructed 
Wc events, which is less than unity due to dilution arising 
from hadronic decays in flight and hadrons misidentificd 
as soft leptons. The terms A and S, which are derived 
from a Monte Carlo (MC) simulation of Wc events and 
the detector response, specify the unfolding from the ob- 
served same-sign subtracted Wc event yield to the mea- 
sured cross section. The cross section is defined through 
A to correspond to the production of a W boson over 
the entire kinematic range associated with a single charm 
quark with pt c > 20 GeV/c, |t7 c | < 1.5. The phase space 
of the charm is restricted to approximately match the de- 
tector acceptance of the charm quark, which minimizes 
the theoretical uncertainties on A. In the determination 
of A, the Wc signal is defined to include events with a 
single charm quark and allows for additional jets; contri- 
butions from all sources of W bosons associated with cc 
pairs are not considered in the acceptance since they can- 
cel out in the same-sign subtraction, owing to the largely 
charge-symmetric detector response. 

The CDF II detector is described in detail elsewhere 
0- The data sample, produced in pp collisions at 
= 1.96 TcV during Run II of the Fermilab Teva- 
tron, was collected between March 2002 and March 2009. 
This analysis is based on an integrated luminosity of 
4.3 ± 0.3 fb" 1 . Events are selected with an inclusive- 
lepton online event selection (trigger) requiring an elec- 
tron (muon) with E T > 18 GeV (p T > 18 GeV/c). The 
trigger and lepton identification efficiencies are measured 
using Z boson decays to electrons and muons, and vary 
between 80% and 97% ||. Further selection requires ex- 
actly one isolated electron (muon), both with isolation 
parameter / < 0.1 [gj], with Et {pt) greater than 20 GeV 
(20 GeV/c) and |?7| < 1.1. The event must also have 
$ T > 25 GeV and exactly one jet with E T > 20 GeV 
and \rj\ < 2.0. The transverse mass of the W boson can- 
didates is required to be greater than 20 GeV/c 2 [Toj | . 
Jets are identified using a fixed-cone algorithm with a 
cone opening of AR = \J (Ary) 2 + (Acf>) 2 = 0.4 and are 
constrained to originate from the pp collision vertex. The 
jet energies are corrected for detector response, multiple 
interactions, and uninstrumented regions of the detector 



Muon candidates inside jets are identified by match- 
ing the trajectories of charged particles (tracks) of the 
jet, as measured in the inner tracking system, with track 
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segments in the muon detectors. An SLT^ |8|, |12j must 
have pt > 3 GeV/c and be within AR < 0.6 of a jet 
axis. Soft electrons from semileptonic heavy-flavor de- 
cay (SLT e ) are identified by tracks with pt > 2 GeV/c 
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that are associated with an electromagnetic shower in the 
central electromagnetic calorimeter, and must lie within 
Ai? < 0.4 of a jet axis. Furthermore, finely segmented 
wire and strip chambers are used to identify the colli- 
mated shower of the electron within the broader hadronic 
shower of the jet. Additional variables to discriminate 
soft electrons are based on the energy deposition, trans- 
verse shower shape, and track-shower distance fl3l [lij]. 
To reduce background from T/7* — ► and Z/j* — > fiji, 
events are discarded where the invariant mass computed 
from the oppositely charged soft muon and primary muon 
is in the 8-11 GeV/c 2 or 70-110 GeV/c 2 ranges. Events 
are also discarded if the jet tagged by a soft muon has an 
electromagnetic fraction greater than 90%, to reduce the 
contamination from Z — > fifj, decays with final-state ra- 
diation off one muon. For the sample with a soft electron 
and a candidate W —> ev decay, the events for which the 
invariant mass between the two electrons is greater than 
45 GeV/c 2 are rejected in order to reduce the background 
from Z/'f* — > ee events. To suppress QCD multijet back- 
ground, we reject events for which the azimuthal angular 
difference between the f, T and the jet is less than 0.3 
rad. In the SLT M (SLT e ) channel we find a total of 1482 
(2494) OS events and 1024 (2088) SS events that satisfy 
the above selection criteria, corresponding to a same-sign 
subtracted excess of 458 (406) events, respectively. 

The dominant backgrounds to Wc are due to the as- 
sociated production of jets with the W boson (W + 
jets, excluding the Wc under investigation), and from 
Drell-Yan production of Z/j*, with and without addi- 
tional jets. Multijet QCD events and small contribu- 
tions from diboson, single top, and tt production are also 
present. Backgrounds are estimated using a combina- 
tion of MC simulation and control regions from the data. 
The MC simulations of W + jets and Zh* + jets pro- 
cesses are performed using ALPGEN (v2.1 [l5|) interfaced 
with PYTHIA (v6.3 [16j) for the parton shower (PS) evo- 
lution. The simulation of the Wc signal is performed 
similarly and is referred to as LO + PS. Modeling of 



heavy-quark hadron decay is provided by evtgen [17 1. 
All samples are simulated using the CTEQ5L PDF sets, 
with Tune BW [3] to model the underlying event and 
the hadronization parameters. Events with a Z — > tt 
decay are also simulated, as well as Zbb and Zee final 
states. The production of Z/j* + jets in the simulation 
is normalized by the measured exclusive Z + 1 jet cross 
section [l9(. 

The W boson events that can mimic the Wc signa- 
ture consist of a W boson associated with heavy-flavor 
quark pairs (bb and cc) or light-flavor (LF) jets. How- 
ever, since this measurement is sensitive to the excess of 
OS over SS events, backgrounds that lead to an equal 
amount of SS and OS events will not contribute to the 
OS excess ascribed to Wc. Therefore, while the recon- 
structed Wbb and Wcc events that pass the selection re- 
quirements are approximately as many as those from the 



Wc signal, they are nearly completely charge symmet- 
ric, since the soft lepton can come from either the b (c) 
or b (c), and the excess of OS over SS events from Wbb 
and Wcc is negligible. On the other hand, W + LF 
events enter the data sample when the jet is identified 
as a charm jet via a misreconstructed soft electron or 
soft muon tag ("mistagging"). Since the same process 
that leads to Wc production describes Wu, a small anti- 
correlation between the charge of the W boson and the 
charge-sign of the tracks in the jets recoiling against the 
W is also observed, leading to a residual background con- 
tribution. We rely on a combination of MC simulations 
and data-driven techniques to estimate this contribution 
to the tagged sample: First the number of W + jets 
events (~97% of which is W + LF) is estimated in the 
sample of events before tagging the jet ( "pretag sample" ) 
by subtracting from the data the initial pretag estimate 
of the signal and all other backgrounds. The number of 
tagged W + jets events is obtained from this pretag es- 
timate using a mistag probability parametrization. The 
procedures for measuring the contribution from real and 
misidentified lepton tags for both the SLT» and SLT e 
were developed in previous data analyses [1, [l3| ; calibra- 
tions were derived using independent samples, validated 
in bb events, and used for the tt cross section measure- 
ments. The probability of misidentifying a hadron as an 
SLT M , denoted as the SLT M mistag probability, is mea- 
sured using a data sample of pions, kaons, and protons 
from A and D* + -tagged D° decays. It is parametrized 
as a function of the track curvature and 77, and is shown 
to describe within ±5% the number of false SLT M tags in 
candidate light-flavor jets of QCD multijet and 7 + jet 
events. A fraction (0.34 ± 0.02)% of the W + LF pretag 
events is expected to be mistagged with a soft muon. The 
SLT e mistag probability is estimated using MC simula- 
tion and checked on a data sample of light-flavor jets from 
QCD multijet events, and a fraction (0.89±0.07)% of the 
W + LF pretag events is expected to be mistagged with 
a soft electron. Both these uncertainties include the sys- 
tematic uncertainty from the corresponding mistag prob- 
ability parametrization. The SLT e mistag probability is 
higher than that of SLT^ also due to contributions from 
photon conversions. The charge asymmetry in W + LF is 
determined directly from the data by applying the mistag 
parametrizations to tracks in the W + jet pretagged sam- 
ple; for the SLT e , the asymmetry is then corrected for 
the 20% contribution to W + LF tags from conversions, 
which is charge symmetric and does not affect the charge 
correlation of tracks in the pretagged sample. The result 
is in agreement with estimates using a simulation of W 
+ LF jets. 

The second largest background to Wc is due to the 
misreconstruction of Z/j* + jets events. The two lep- 
tons from the Z/^f* decay can be misidentified as one 
lepton from a W boson decay and one soft lepton, re- 
sulting in approximately 90% charge asymmetry. These 
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events are suppressed by the veto on the if-mass region. 
Alternatively, only one lepton from the Z boson decay 
is reconstructed in the event, which is typically assigned 
to be a VF-decay lepton. In this case, the soft lepton 
results from the decay of heavy flavor or from the mis- 
reconstruction of a track from hadrons, and these events 
carry approximately 40% asymmetry. The overall aver- 
age charge asymmetry of Z/j* + jets for SLT e is smaller 
than for SLT^ because of the stricter requirements on the 
dielectron mass. 

Events due to QCD multijet production can enter the 
selection through hadronic misidentification or heavy- 
flavor decay. Missing transverse energy can arise from 
mismeasured jet energy, detector effects, or neutrinos in 
the decay chain. We estimate this background by releas- 
ing the missing energy requirement on the events and 
fitting templates of the $ T distribution for the QCD mul- 
tijet component, separately for OS and SS events. The 
template distribution for QCD multijet events is derived 
from a jet-enriched data sample in which candidate elec- 
trons fail two of the electron identification criteria. The 
remaining sample composition is modeled with MC sim- 
ulations. 

Finally, the production of dibosons (WW, WZ, ZZ) 
and tt is modeled with a pythia (v6.4) MC calcula- 
tion, while single top-quark production is simulated using 
MADEVENT [20j. These processes together contribute a 
few percent of the total background, and their production 
cross sections are well established. The WW events con- 
tribute the most and have a strong charge asymmetry. 
Table [J summarizes the data and the estimated back- 
ground. 

We assume that the total OS-SS rates observed in 
the data, after subtracting the background contributions 
listed in Table U are due to the Wc signal; the SS- 
subtracted rates for the signal are then 287 ± 50(stat) ± 
32(syst) and 149 ± 68(stat) ± 26(syst) events, for the 
SLT M and SLT e tagged samples, respectively. This is 
higher than the number of OS-SS events predicted by 
the LO + PS simulation of Wc production. Figure Q] 
shows the distributions of the measured px spectrum 
for SLT muons and electrons in tagged events, com- 
pared to the predicted spectrum given by the Wc signal 
and estimated background. For each contribution, SS 
events are subtracted. The Wc production cross section 
is calculated using Eq. ([1]), with awe = Cw/+ E + &w-c-> 
B{W -> iv) = 0.108 ± 0.009 [2l|, pre > 20 GeV/c, and 
\r) c \ < 1.5; the values of the dilution S for Wc events are 
given in Table [I] We measure owe X B (W — > iv) = 13.4 ± 
2.3(stat)±|^(syst)±^(lum) pb axida Wc xB(W -» iv) = 
15.0 ± 6.8(stat)t£g(syst) ± 1.2 (lum) pb from the SLT M 
and SLT e samples, respectively. 

Systematic uncertainties are shown in Table [TTJ The 
uncertainty on the SLT tagging includes contributions 
from the measurements of the efficiency of tagging lep- 
tons in a jet environment and of mistagging [8j, [l3|. The 



TABLE I: Summary of data and backgrounds in the SLT M - 
tagged and SLT e -tagged W + 1 jet samples. The expected 
Wc contribution is shown from a LO + PS calculation. 



Source 



Events Asymmetry OS-SS 



SLT. 



W + LF, bb, cc 


1808 ±271 


0.05 ±0.01 


86 ± 14 


Z/l* + jets 


132 ±30 


0.63 ±0.02 


84 ± 18 


QCD multij. 


308 ± 17 


-0.03 ±0.07 


-8 ±17 


Diboson, t(i) 


26 ±3 


0.33 ±0.01 


9± 1 


Wc (LO + PS) 


214 ±19 


0.75 ±0.03 


161 ± 13 


Total expected 


2488 ± 274 




331 ±37 


Data 


2506 




458 


SLT e 


W + LF, bb, cc 


4076 ± 305 


0.04 ±0.01 


174 ± 19 


Z/l* + jets 


138 ±29 


0.26 ±0.01 


36 ±7 


QCD multij. 


374 ± 12 


0.07 ±0.03 


27 ± 12 


Diboson, t(t) 


35 ±3 


0.58 ±0.01 


20 ±2 


Wc (LO + PS) 


174 ±16 


0.45 ±0.02 


78 ±7 


Total expected 


4797 ± 307 




336 ± 28 


Data 


4582 




406 



TABLE II: Summary of systematic uncertainties, as a per- 
centage of the measured Wc cross section. Numbers shown 
in bold font indicate uncertainties treated as uncorrelated in 
the combination of the channels. 



Source 


SLT^ 


SLT e 


SLT uncertainties 


±9.2 


±16.6 


QCD multijet estimate 


±6.3 


±9.9 


Initial and final state radiation 


±6.0 


±6.0 


Background cross sections 


±5.7 


±4.7 


c quark hadronization 


±4.6 


±4.6 


PDFs 


±3.6 


±3.6 


VL-lepton ID 


±2.2 


±2.2 


Jet energy calibration 


±2.0 


±2.0 


Factorization, renormalization scales 


±1.3 


±1.3 


Total 


±15.4 


±21.8 


Luminosity 


±7.9 


±8.3 



uncertainty on the backgrounds includes contributions 
from the theoretical cross sections, from the estimation 
technique, and from statistics for the backgrounds evalu- 
ated with inputs from a data control region. For the Z/"f* 
background, the dominant uncertainty on the event yield 
estimate comes from the measured Z cross section uncer- 
tainty. To measure the effects of initial- and final-state 
gluon radiation, we measure the Wc acceptance in differ- 
ent sam ples with the radiation enhanced or reduced, as 
in Ref. [22j . We compare charm jets modelled with the 
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FIG. 1: The soft muon and soft electron pr distributions. 
The Wc contribution is derived from a LO + PS simulation 
normalized to the measured cross section. 



pythia and herwig [23|, l2J] MC calculations to evalu- 
ate the uncertainty due to different hadronization mod- 
els. The PDF uncertainty is derived by remeasuring the 
acceptance using the CTEQ and MRST sets, following 
the same prescription as in Ref. |22j. The MC model- 
ing of the efficiency for identifying the leptons from the 
W boson decay ("W lepton ID") is measured using Z 
boson data and MC samples. The charge misidentifica- 
tion rate is less than 1% and therefore has a negligible 
effect. The uncertainty due to the jet energy calibration 
(JES) is measured by shifting the energies of the jets in 



111. The 



the Wc MC simulation by ±ler of the JES 
certainty on the acceptance due to the factorization and 
renormalization scales is estimated by varying them in 
the ALPGEN MC between 1/2 and twice the transverse 
mass of the W boson, as well as using the charm quark 

PT- 

The results from the two SLT-tagged samples are 
combined by performing a profile likelihood ratio 
minimization [25]. Systematic uncertainties are in- 



cluded as nuisance parameters with Gaussian con- 
straints whose widths are fixed to the respective un- 
certainties, and are assumed to be either fully cor- 
related, if they are shared between the two chan- 
nels, or uncorrelated if not. The combination of the 
cross section yields ay/c (pre > 20 GeV/c, \r] c \ < 1.5) x 
B{W -> iv) = 13.6±2.2(stat)t 2 j(syst)±l.l(lum) pb = 
13.6^3 * pb. The significance for the Wc signal is derived 
from the ratio of profile-likelihoods A, with — 21nA in the 
hypothesis of no signal being interpreted as following a 
X 2 -distribution, and is calculated to be 5.7cr. The mea- 
surement is in agreement with a NLO calculation over 
the same phase space of 11.4 ± 1.3 pb [26[, where the 
renormalization and factorization scales have been set to 
half the W boson mass, and varied between 5 GeV and 
80 GeV in the uncertainty. The uncertainty also includes 
PDF variations using the CTEQ6M and MSTW2008 
sets. The result can be also compared to the LO pre- 
diction of 8.2 ± 1.5 pb (2(|, giving a measurement to LO 
cross section ratio for this kinematic region of 1.6 ± 0.5. 
Since the majority of Wc production proceeds through c 
to s quark coupling, we can relate the measured value of 
the cross section with the thoretical prediction and de- 
rive \V CS \. Using a% e c orv = 9.8(±1.1)|U CS | 2 + 2.1(±0.2) pb 
[26| we obtain |T4 S | = 1.08±0.16, where the uncertainties 
in the cross section measurement and in the theoretical 
prediction have been added in quadrature. Restricting 
the range of \V CS \ to the interval [0,1], a lower limit of 
\V CS \ > 0.71 at the 95% confidence level is extracted. 

In conclusion, we present the first observation of pro- 
duction of a W boson associated to a single charm 
quark, with a significance of 5.7 standard deviations. 
This is obtained in data corresponding to 4.3 fb _1 pro- 
duced in pp collisions at y/s = 1.96 TeV and collected 
by the CDF experiment. The measured cross section 
is ex Wc (pre > 20 GeV/c, |r? c | < 1.5) x B (W -» iv) = 
13.6^3 i pb, which is in agreement with a SM NLO cal- 
culation. A direct determination of |V CS | = 1.08 ± 0.16 
and the 95% C.L. lower limit of V cs \ > 0.71 are extracted 
assuming that the Wc production through c to s quark 
coupling is dominant. 
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